Quantifying, characterizing and controlling information flow in ultracold atomic gases 
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We study quantum information flow in a model comprising of an impurity qubit immersed in a 
Bose-Einstein condensed reservoir. We demonstrate how information flux between the qubit and the 
condensate can be manipulated by engineering the ultracold reservoir within experimentally realistic 
limits. We place a particular emphasis on non-Markovian dynamics, characterized by a reversed flow 
of information from the background gas to the qubit and identify a controllable crossover between 
Markovian and non-Markovian dynamics in the parameter space of the model. 
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Introduction- In the last decades high precision con- 
trol of ultracold atomic gases has allowed the realization 
of beautiful experiments unveiling fundamental phenom- 
ena in the physics of many-body quantum systems at 
low temperatures. Key examples are the observation of 
Anderson localization the superfluid-Mott insulator 
transition Q , the creation of Tonks-Girardeau gases [j| , 
and the atom laser [ij], just to mention a few. 
More recently, hybrid systems composed of quantum 
dots, single trapped ions and optical lattices coupled 
to Bose-Einstein condensates (BECs) have been studied 
both theoretically and experimentally [j| . These systems 
are studied in the framework of open quantum systems 
0, effectively described as one or more two-level sys- 
tems (qubits) interacting with a reservoir consisting of 
the ultracold gas. The possibility of manipulating cru- 
cial parameters of the reservoir, such as the scattering 
length 7], combined with the continuous improvements 
in quantum control of qubits highlights the enormous po- 
tential of hybrid systems as quantum simulators of both 
condensed matter models and open quantum systems. 
Lately, the open systems community has given substan- 
tial interest to the dynamics of quantum information flow 
due to several proposals to link it to the division of quan- 
tum processes into Markovian and non-Markovian ones 
8 111. The latter ones have been defined as processes 
where an open system recovers some previously lost infor- 
mation and therefore temporarily combats the destruc- 
tive effect of the environment as a sink for quantum prop- 
erties. Such processes are very desirable from a quantum 
information processing point of view, since they increase 
the operational time of qubits [121 ] . Sophisticated con- 
trol over detrimental environmental effects enabled by, 
for example, hybrid quantum systems will be very bene- 
ficial for state-of-the-art quantum processors. 
Moreover, understanding of the decoherence processes in 
hybrid systems gives insight into the basic physical mech- 



anisms ruling the interaction of the qubit and the cold 
atomic cloud surrounding it. It is also a key pre-requisite 
to the realization of quantum simulators since it allows 
for the identification of the relevant parameters for on- 
demand control of the exchange of quantum information 
between a qubit and its ultracold environment. 
We present here the first characterization of non- 
Markovian effects in the context of ultracold gases. More 
specifically, inspired by a model introduced in Ref. fl3j j . 
we consider an impurity atom in a double well potential, 
i.e., an effective qubit system, interacting with a BEC 
environment. We study information flow using the non- 
Markovianity quantifier presented in Ref. [|[ and show 
how information flux can be manipulated by experimen- 
tally achievable means such as changing the scattering 
length, the effective dimension of the background gas or 
the trapping geometry of the qubit. We find a crossover 
between Markovian and non-Markovian dynamics in the 
parameter space of the model and uncover the physical 
mechanisms at the root of non-Markovian phenomena 
induced by the ultracold background gas. Our findings 
pave the way to the realization of quantum simulators 
for non-Markovian open quantum system models with 
ultracold atomic gases. 

The model- We consider an impurity atom trapped in 
a deep double well potential Va(i"). The impurity atom 
forms a qubit system with the two qubit states repre- 
sented by the occupation of the impurity atom in the left 
or the right well, \L) and \R) respectively. The impurity 
atom couples to a bosonic background gas B trapped in 
a shallow potential Vs(r), which forms a Bose-Einstein 
condensed environment for the qubit system. The Hamil- 
tonian for this system, derived in Ref. is 
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where a, = \R)(R\-\L)(L\ and E k = ^ e k [e k + 2g ( ° ] "no] 
is the energy of £;-th Bogoliubov mode Ck of the con- 
densate with boson-boson coupling frequency g^ and 
condensate density nn- D denotes the effective dimen- 
sion of the environment. The energy of a free mode is 
£k = h 2 k 2 1 (2ms) where k = |k| and tub is the mass of 
a background gas particle. Furthermore, <7k and £k are 
coupling constants that depend on the spatial form of the 
states \L) and \R) and on the shape of the Bogoliubov 
modes. Their specific form is elaborated in Ref. (l3| . 
When the background gas is at zero temperature the re- 
duced dynamics of the impurity atom is captured by the 
following time- local master equation (ME): 



dp(t) 
dt 



= A(t)[a z ,p]+ 1 (t)[a zP (t)a z - \{u z a z , p{t)}]. (2) 



Quantity A(i) renormalizes the energy of the qubit but 
has no qualitative effect on the dissipative dynamics. In- 
stead in this work we are interested in the decay rate 
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where gAB is the impurity-boson coupling frequency, r is 
a trap parameter, and L is half the distance between the 
two wells of the double well potential. 
We have derived ME © using the time-convolutionlcss 
projection operator technique to second order in the 
coupling constant gAB 6]. Remarkably, in this case 
the second order ME describes the reduced dynamics 
exactly [14| . Solving the ME reveals that the impu- 
rity atom dephases without exchanging energy with the 
background gas. More precisely, pu(t) = p„(0) and 
Pij{t) = e~ r Wpy(0) when i ^ j, where p i} = (i\p\j), 
i,j = R,L. The decoherence function T(t) — f*dsj(s) 



coincides with that derived in Ref. 13[ , however here we 
wish to stress the connection between the decay rate and 
the non-Markovian features. The authors of Ref. [3| dis- 
covered situations when the decoherence function T(t) is 
non-monotonic and conjectured that this is due to non- 
Markovian effects in the reduced dynamics. Already the 
form of the ME ([2]) supports this intuition; the theory of 
non-Markovian quantum jumps has shown that there is a 
profound connection between non-Markovian effects and 
negative regions of the deca y ra tes of Lindblad-structured 
MEs as the one of Eq. ^ [l5| . In the following we con- 
firm that this is indeed true for this model and, moreover, 
expose the physical mechanisms at the root of this non- 
Markovian phenomena. 

Non-Markovianity measure- Breuer, Lainc and Pi- 
ilo (BLP) have proposed a rigorous definition for non- 
Markovianity of a quantum channel $ based on the 
dynamics of the so-called information flux cr(t) = 
dD[pi(t), p2(t)]/dt 0. This is the temporal change in 



the distinguishability D[pi(t) , p 2 (t)] = — y»2 (*) 1 1 1 

of two evolving quantum states pi, 2 (t) = §(t)pi t2 (0) as 
measured by the trace distance. Negative information 
flux describes information leaking from the system to its 
environment and it is associated to Markovian dynamics. 
Instead, if it is possible to find a pair of states pi, 2 {0) 
for which the information flux is positive for some in- 
terval of time, that is, the system regains some of the 
previously lost information, then process <£> is consid- 
ered non-Markovian. The amount of non-Markovianity 
is defined to be the maximal amount of information that 
the system may recover from its environment, formally 
Nblp = max pi 2 J a>Q dsa(s). 

For the model studied in this Letter we find that a(t) > 
if and only if 7(i) < 0, that is, the process is non- 
Markovian precisely when the decay rate can take tem- 
porarily negative values. Within experimentally relevant 
values of the physical parameters we have discovered at 
most a single time- interval t € [a,b] when the decay rate 
is negative and information flows back to the system after 
an initial period of information loss. Therefore, instead 
of using the original measure faithfully and quantifying 
non-Markovianity as the maximal amount of information 
that the system may recover, we introduce a normalized 
quantity that reveals the maximal fraction of the previ- 
ously lost information that the system can recover: 
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Unlike Nblp, the modified quantifier J\f is bounded be- 
tween zero (system only leaks information) and one (sys- 
tem regains all previously lost information) and is there- 
fore more meaningful as a number. We have confirmed 
numerically that in the relevant case of dephasing noise 
the above quantity is maximized for the same pair of 
initial states that maximize Nblp- These are the states 
whose Bloch vectors lie on the opposite sides of the equa- 
tor of the Bloch sphere [l6|. Using these states in the 
general expression of Eq. (U) we find the analytic expres- 
sion of the non-Markovianity measure for a dephasing 
qubit to be 

e -r(fc) _ e -r(a) rt 
Nd -P h = e -r(o) _ e -i» ' r W = J o ds ^- ( 5 ) 

We are now ready to study how changes in the back- 
ground scattering length and in the dimensionality of the 
BEC affect the dynamics of information flow. 

Three-dimensional BEC- As a first step we consider 
a 3D background BEC with equal confinement of the 
background gas in all directions. We consider a 87 Rb- 
condensate of density n$ = tiq = 10 20 m -3 and 23 Na im- 
purity atoms trapped in an optical lattice with lattice 
wavelength A — 600nm and trap parameter r = 45nm. 
The impurity-boson coupling is gAB = ^ttH 2 ciab /m-AB, 
where tuab = niAiriB / '(rriA + rns) and itia and tub are 
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FIG. 1. (Color online) Non-Markovianity measure J\fdeph as 
a function of the scattering length of the background gas as 
with values of well separation L — 50nm (red dashed line), 
L = 75nm (blue dotted line) and L — lOOnm (black solid 
line). The inset shows Mde. v h for very small values of the 
scattering lenght. 



the masses of the impurity atoms and the bosons, re- 
spectively, and clab = 55 ao, where ao is the Bohr ra- 
dius. Similarly the boson-boson coupling frequency is 
= AtxT) 2 <xb Ivtib but now we assume that the s-wave 
scattering length of the background gas can be tuned 
from its natural value as = a-m ~ 5.3nm via Feshbach 
resonances. We explore a range of values of as consistent 
with the assumption of dilute gas and with the regime of 
weakly interacting gases. The latter is a stronger condi- 



tion, requiring y/ a^rio < 1. As a consequence, we can 
tune the scattering length up to a maximum value given 
by a B « 3a Rb . 

Figure. [1] shows the non-Markovianity measure Afdeph as 
a function of as for three different values of the well 
separation L. Increasing L magnifies the fraction of 
recovered information flow due to the increased ability 
of the condensate to resolve the qubit system. Simi- 
larly, non-Markovian effects are amplified for stronger 
interaction of the condensate. However we find that 
the scattering length alone plays a crucial role in the 
emergence of non-Markovian reservoir memory effects. 
When the background gas is free or very weakly inter- 
acting, < as < a c £ %t ~ 0.034 a^, the qubit only 
leaks information to the BEC environment. Instead for 
a strong enough interaction strength of the background 
gas, as > af', the condensate can take on the role of 
information storage and feed some information back to 
the qubit. This result holds for any value of L. This 
finding challenges the conclusion of Ref . [l3[ , where the 
scattering length dependent Markovian/non-Markovian 
crossover was only attributed to the ID case. We have 
discovered that the situation is indeed more subtle and 
we will show next that the crossover point exists in all 
three dimensions. 



Lower dimensions- By a suitable modification of the 
potential of the condensate Vb (r) we can create a quasi- 
2D background gas where the gas is trapped in a slightly 
anisotropic, pancake-shaped harmonic trap. Assuming 
that the scattering length is still much smaller than the 
axial length of the condensate, as -C a z , the coupling 
term is modified to g 2 ^? = v / 87ra 2 aFj/(m fja ? ) and the 
2D condensate density is n% = \f~iinQa z [171 ]. Within 
the limits of a dilute gas we can increase the scatter- 
ing length up to as ps 2auh- The potential Vb{t) can 
be also modified to create a cigar-shaped quasi-lD back- 
ground gas with transversal width a±. The consequent 
coupling is — 2Ti 2 clb / (rnB a \) and the ID density is 
n i = n o 7ra j_i again provided that gas is weakly enough 
confined, ag -C aj_ [18]. In the quasi-lD regime dilutc- 
ness of the gas allows at most as ^ am- 
In Fig. [2] we plot the non-Markovianity measure Ndeph 
in the quasi-lD, quasi-2D and 3D cases. In lower di- 



mensions we find the critical values 



crit, 2D 



0.122 ai » 



and 



crit, ID 
B 



0.183a/». Clearly when the dimensional- 



ity of the background gas is lowered the cross-over value 
of the scattering length a^ lt increases. Crucially, as we 
remarked before, a^ l *' 3D > and therefore it is possible 
to create both Markovian and non-Markovian dynamical 
processes in all three dimensions. 

We note here that the quantities we have chosen to vary, 
namely, the scattering length a b and the well separation 
L, are indeed the most relevant quantities for manipu- 
lating the information flow-back. The trap parameter 
t determines the trapping frequency of the double well 
trap Vk(r) and acts as a natural cut-off parameter in 
the decay rate of Eq. Q. As long as the double well 
trap is deep enough to prevent hopping between the two 
sites the particular value of r has only a minor effect on 
information flow. It is also clear from the form of the 
decay rate that the boson-impurity coupling §ab cannot 
affect the Markovian/non-Markovian crossover. More- 
over, we have found that its value has negligible effect 
on the non-Markovianity quantifier Ndeph ■ Intuitively, it 
affects the amount of outgoing and incoming information 
almost equally but leaves their ratio unchanged. In order 
to explain the key non-Markovian features in the dynam- 
ics of the qubit system when the dimensionality and the 
scattering length of the background gas vary, we need to 
take a closer look at the spectrum of the BEC reservoir. 

Environmental spectrum- The crossover between 
Markovian and non-Markovian processes is best under- 
stood in terms of the environmental spectrum J(u>). 
Consider the general dephasing model introduced by 
Palma, Suominen and Ekert describing qubit dynam- 
ics: pu(t) = Pii(0) and pij(t) = e" r(t) p y (0), where 
f(t) ~ J dojJ(uj)(l — cos tot) /lu The dynamical pro- 
cess is non-Markovian if and only if j(t) = dT(t)/dt < 
for some interval of time. We assume an Ohmic-type 
spectrum J(u>) ~ lo s and recall the convention that the 
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FIG. 2. (Color online) Non-Markovianity measure J\fdeph 
a function of the scattering length of the background gas as 
when the background gas is three dimensional (red dashed 
line), quasi-two dimensional (blue dotted line) and quasi-one 
dimensional (black solid line) . The inset shows a longer range 
of the scattering length as- In all figures the well separation 
is L = 75nm. 



spectrum is sub-Ohmic when s < 1, Ohmic when s = 1 
or super-Ohmic when s > 1. Introducing an ad hoc ex- 
ponential cut-off so that J{lu) — uj s exp{— ui 2 /uiq}, where 
loc is the cut-off frequency, it is straightforward to show 
that the dynamics is non-Markovian when s > s crtt — 2. 
Therefore, in a general setting, a qubit dephasing under 
the effect of either a sub-Ohmic or an Ohmic environ- 
ment can only leak information to its environment. If 
the environment has a super-Ohmic spectrum the issue 
is less straightforward: only if the spectrum is sufficiently 
super-Ohmic with s > s crlt , information can flow back 
to the system from the environment. 
The qubit in an ultracold bosonic environment consid- 
ered in this work is a special case of the model above 
with J(uj) = 2k lffk| 2 (5(^— Ek). The reservoir spectrum 
J (to) is very complex due to the complicated form of the 
coupling constant g^. However, it can be shown that in 
the case of a free background gas in one, two or three 
dimensions the spectrum is sub-Ohmic, Ohmic or super- 
Ohmic, respectively [3]. The spectrum changes criti- 
cally when one considers the boson-boson coupling quan- 
tified by the scattering length a^. In this case increasing 
the scattering length effectively increases the value of s. 
Hence when we increase as in the ID case, the spectrum 
changes from sub-Ohmic to Ohmic to super-Ohmic and 
once a critical threshold of super-Ohmicity is reached the 
environment can feed information back to the system. In 
the 2D non-interacting case the spectrum is Ohmic and 
a weaker interaction is required to reach the crossover 
point s crlt , leading to a C g lt ' 2B < a C g lt ' 1D . Finally, in the 
3D case the spectrum is already super-Ohmic in the non- 
interacting case, although not super-Ohmic enough to 
give rise to non-Markovian dynamics. Already a small 
increase in the scattering length modifies the spectrum so 



that the direction of information flow can be temporarily 
reversed. 

Conclusion- We have studied quantum information flux 
in an ultracold hybrid system of an impurity atom im- 
mersed in a BEC environment. We have shown explicitly 
how precise control of the ultracold background gas af- 
fects the spectrum felt by the qubit and therefore enables 
the manipulation of the qubit dynamics and the informa- 
tion flux. In particular, we have discovered experimen- 
tally accessible means to reach non-Markovian dynamical 
regimes, where the background gas may feed informa- 
tion back to the qubit instead of acting only as a sink 
for information. Such quantum reservoir engineering is 
fundamental for understanding decoherence processes in 
quantum information processing and, more specifically, 
for the realization of quantum simulators. 
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